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Abstract 

We present the first calculation of the two-loop electroweak fermionic correction 
to the flavour-dependent effective weak-mixing angle for bottom quarks, sin^ 9^. 
For the evaluation of the missing two-loop vertex diagrams, two methods are em- 
ployed, one based on a semi-numerical Bernstein-Tkachov algorithm and the second 
on asymptotic expansions in the large top-quark mass. A third method based on 
dispersion relations is used for checking the basic loop integrals. We find that for 
small Higgs-boson mass values, Mh oc 100 GeV, the correction is sizable, of order 
0(10-4). 
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1 Introduction 



Experiments at LEP, SLC and Tevatron have provided a large number of high-precision 
data, which, being supplemented by detailed studies of higher-order corrections, allow to 
probe the Standard Model at the loop level and subsequently to predict the mass of the 
Higgs boson. In this context, the leptonic effective weak-mixing angle, sin^ ^efF plays the 
most crucial role. It can be defined through the effective vector and axial- vector couplings, 
gy and g''^, of the Z boson to leptons (/) at the Z-boson pole. 
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Singer = 7 1 + (1) 



The effective weak-mixing angle can be related to the on-shell Weinberg angle, sin^ 9^, as 

sin^^if = sin2e^«:, (2) 

where sin^ 6^ = 1 — M^/M| and k = 1+Ak. At tree level. An = and sin^ = sin^ 9w 
The form factor An incorporates the higher-order loop corrections. Usually, the Vl^-boson 
mass, Mw, is not treated as an input parameter but it is calculated from the Fermi 
constant, G^, which is precisely known from the muon lifetime. The relation between 
Mw and can be cast in the form 

where the quantity Ar [1] contains all higher-order corrections. The presently most ac- 
curate calculation of the ly-boson mass includes full two-loop and leading higher-order 
corrections [2]. On the other hand, the quantity k in Eq. ([2]) incorporates all corrections 
to the form factors of the Zll vertex. Recently, the calculation of the two-loop electroweak 
corrections has been completed [3-7]. The uncertainty on sin^ 6^^^^ due to unknown higher 
orders has been estimated to be 0.000047, which is substantially smaller than the error of 
the current experimental value sin^ 9^^^^ = 0.23153 ± 0.00016 [8], but still larger than the 
expected precision, 1.3 x 10~^, of a future high-luminosity linear collider running at the 
Z-boson pole [9]. 

The experimental value for sin^ 9^^^^ is determined from six asymmetry measurements, 
A^PQ, Ai{Pt), ^;(SLD), A^pq, A^jf^, and Q^^- Of those, the average leptonic and hadronic 
measurements differ by 3.2 standard deviations, which is one of the largest discrepancies 
within the Standard Model. The main impact stems from two measurements, the left- 
right asymmetry with a polarised electron beam at SLD, /j, and the forward-backward 
asymmetry for bottom quarks at LEP, ^^g. On the experimental side, the only possible 
source of this discrepancy are uncertainties in external input parameters, in particular 
parameters describing the production and decay of heavy-flavoured hadrons; see Section 
5 of Ref. [8] for a discussion. However, the interpretation of the asymmetry measurements 
in terms of sin^ 9^^^^ requires also some theoretical input. The leptonic asymmetries depend 



2 



on lepton couplings only and can be translated straightforwardly into the leptonic effective 
weak-mixing angle, with small corrections due to s- and t-channel photon exchange. By 
contrast, the hadronic observables, .4^^, A^p% and Q^g, depend on the quark couplings, 
gy^. These couplings are associated with a flavour-dependent hadronic effective weak- 
mixing angle, sin^ 9^^, 

sm'eil = ^^(^l + Re^y q = d,u,s,c,b. (4) 

The forward-backward pole asymmetry of a quark g, A^jf^, is related to the effective 
couplings, gy and g-^, and the effective weak-mixing angle, sin^^^|, b}|3 

3 

ApQ = -AeAq, (5) 

with 

' ig'v)' + {9{? 1 - A\Q;\ sin^ + 8Q2 sin^ qQ' 

At tree level, sin^ 6'^ and sin^ O^^^' are identical, but the relations between these quantities 
receive sizable radiative corrections that need to be included in the analysis. Note that, 
due to the small electric charge of the bottom quark, Qb = —1/3, the parameter At is 
close to 1, and A^J^^ is only weakly sensitive to sin^ 9^. Therefore, it seems unlikely that 
the discrepancy between ^°_r -^fb could be explained by radiative corrections. Nev- 
ertheless, the theoretical prediction for sin^ 9^^ enters in the Standard- Model fits through 
several observables, so that a precise prediction of this quantity is important for a robust 
analysis. 

For all fermions except bottom quarks, the known radiative corrections to sin^ 9^^ in- 
clude at least two-loop fermionic electroweak contributions and some leading higher-order 
corrections; see Ref. [7] for details. However, for the Zbb vertex only one-loop correc- 
tions, leading two-loop corrections for large values of the top-quark mass of 0{a'^m^), 
and two- and three-loop QCD corrections have been calculated [10] and included in the 
Zfitter program [11] (see also the new program Gfitter [12]), which is widely used for 
global Standard-Model fits. The remaining two-loop electroweak corrections beyond the 
0{a'^m^) contributions are still unknown, although they are expected to be larger than 
the 0{a^mf) term, based on experience from sin^ ^efF *• As a result, the present treatment 
of higher-order electroweak corrections leads to inconsistencies, for example in since 
the corrections to sin^ and Ae include two-loop and leading three-loop corrections 
that are absent for sin^ 6'^^ and Ab (see recent discussion in Ref. [13]). 

In this paper, the part of the missing two-loop corrections to sin^ 9^^ with closed 
fermion loops is presented. We begin by explaining the techniques employed for the 



^ Owing to the non-zero bottom-quark mass, the Zbb vertex also has a scalar part, besides the vector 



and axial-vector parts. We checked explicitly that the contribution of this scalar form factor to A'^p'g is 
more than a factor 1000 smaller than the current experimental uncertainty and thus truly negligible. 
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calculation in the next section. In Section [3l numerical results for sin^ 9^^ are given before 
the summary in Section HI 

2 Outline of the calculation 
2.1 General approach 

We work in the Standard Model and adopt the on-shell renormalisation scheme, which 
relates the renormalised masses and couplings to physical observables. Details on the 
renormalisation scheme and explicit expressions for the relevant counterterms can be 
found in Refs. [7, 14] . For the loop integrations, we employ dimensional regularisation. 
The problem of 75 matrices in two-loop vertex diagrams with fermion triangle sub-loops is 
treated in the same way as in Refs. [3,4,7], by evaluating the finite non-anticommutative 
contribution from 75 to the vertex diagrams in four dimensions. Most aspects connected 
with the calculation of the effective weak-mixing angle for the Zbb vertex are the same as 
for the leptonic effective weak-mixing angle and are discussed in detail in Ref. [7]. 

The contributions for the two-loop renormalisation terms are identical to the case of 
sin^O^^^^, with the exception of the two- loop bottom-quark wave-function counterterm, 
which involves new self-energy diagrams with internal top-quark propagators; the first 
terms of this quantity are given in Ref. [15]. For the two- loop Zbb vertex corrections, on 
the other hand, a number of new three-point diagrams need to be computed. In general, 
electroweak two-loop corrections can be divided into two groups, which are separately 
finite and gauge invariant: fermionic corrections (with at least one closed fermion loop) 
and bosonic corrections (without any closed fermion loops). In this article, we focus on 
the fermionic diagrams as a first step. 

For the purpose of this calculation, all light-quark masses are neglected in the two- 
loop diagrams, including the bottom-quark mass. As a result, for many diagrams, known 
results from the sin^6'g'j^* calculation can be used [3,4,7]. The loop integrals for dia- 
grams with closed massless-fermion loops are given in analytical form, while large-mass 
expansions were employed for diagrams with top quarks in the loops. 

However, the two-loop corrections to sin^ 9^^ include a new group of integrals that were 
not covered in previous calculations of sin^6'g^*, stemming from diagrams with internal 
ly-boson and top-quark propagators; see Fig. [H The computation of these diagrams will 
be discussed in detail in the following subsections. The two-loop diagrams are computed 
with several independent methods, so that cross checks can be performed. The first 
method, based on the observation that all new diagrams in Fig. [1] include internal top- 
quark propagators, uses asymptotic expansions for large top-quark mass. This method 
was already employed successfully for the calculation of siia^ 9^^^^ [7]. For references on 
the subject, we refer the reader to Ref. [16]. 

Secondly, we develop a code for the evaluation of Feynman diagrams with a semi- 
numerical method, based on the Bernstein-Tkachov (BT) method of Ref. [17]. This 
method had already been used previously for one-loop problems [18]. In a recent series of 



4 



(A) (B) 




Figure 1: Set of Feynman diagrams required for the calculation of the fermionic two- loop 
corrections to the Zbb vertex, but absent in the sin^ 9^^^^ case. Thick solid lines denote 
top-quark propagators, while thin lines represent light fermions. 

papers [19], it was extended to general two-loop vertices, and some applications to two- 
loop problems are already known: the leptonic effective weak-mixing angle was presented 
in Ref . [5] and corrections to the — 77 decay width in Ref . [20] . 

Finally, we use another semi- numerical method based on dispersion relations [21], 
which was also used previously for sin^ 6*1.^* [7]. This method allows us to evaluate all 
self-energy diagrams, the vertex diagrams in Figs. [I](A)-(D), as well as the scalar integrals 
with the topology of Figs. [1](E)-(G). However, due to problems with the complex tensor 
structure, the complete diagrams in Figs.[T](E)-(G) cannot be checked with this technique. 

In the next subsections, we explain the applications of these methods for our purposes 
and present a comparison between them. 

2.2 Asymptotic expansions 

We perform an expansion in a parameter x, where 

Ml 1 

4 ^ ^ 

For any two-loop problem, there are four regions to consider. Let ki and ^2 represent 
the internal momenta in the loops and p stand for any external momentum, while m 
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generically denotes all masses that are small compared to rrit, m < rrit. In our case, 
m = M\Y, Mz- Then the four regions can be identified as follows: 

1) ki ^ rrit and /c2 ~ (expansions in small parameters: p and m) 

2) fci ~ m and k2 ~ rrit (expansions in small parameters: p, ki and m) 

3) ki ~ rrit and k2 ^ m (expansions in small parameters: p, k2 and m) 

4) A;i ~ m and k2 ^ m (expansions in small parameters: p, ki, /c2 and m) 

This method allows us to represent two-loop vertex diagrams by a sum of simpler in- 
tegrals, namely two-loop propagator and vacuum integrals, plus one-loop integrals. How- 
ever, higher orders in the expansion lead to higher powers of propagator denominators in 
these integrals. This is not a problem for one-loop or vacuum integrals, as analytic rela- 
tions are well known; for relations and references, see, for example, Ref. [16]. For two- loop 
propagator integrals, we employ the Laporta algorithm, as proposed in Ref. [22]. This 
algorithm allows us to automatically reduce complicated multi-loop integrals with non- 
trivial numerators to a smaller set of master integrals with unit numerators. In addition 
to the well-known integration by parts relations [23], Lorentz identities [24] can be used 
for faster performance. In our approach, the Laporta integral reduction is accomplished 
with the help of the program IdSolver [4,25]. 

As already observed for the two-loop vertex diagrams in the leptonic case [7], this 
expansion has a fast convergence behaviour. After performing the expansion down to 
fifth order, 0{x^), a precision of 10~^ in the final result for the two-loop part of sin^ 6'^^ 
is obtained. 

2.3 Semi-numerical integration based on the BT algorithm 

Any Feynman diagram can be described by an integral 



where Q(x) and Vi(x) are polynomials of x, Q is the numerator of a Feynman integral, 
Vi is a denominator of the Feynman integral with the power /ij, which depends on e, 
X = (xi, ...,Xn) represents the n-dimensional vector space of Feynman parameters, and S 
is the integration region defined by 



Tkachov proved the existence of an algorithm [17], later called BT algorithm, that can 
transform a Feynman integral of the type in Eq. ([8]) into a form with better arrangement 
of divergences, which is consequently more suitable for numerical integration. However, 
until now, compact-form solutions of this algorithm are only known at the one-loop order, 
which will be described in the following. 

For one-loop cases, V(x) is a quadratic polynomial of x of the form 






(9) 




(10) 
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where W is a, n x n matrix, R is a n-dimensional vector and Z is a scalar number. Then 
one can show that the following relation is fulfilled: 

^ ^ / (x + A)a \ 

where A = (Z-R^W~^R) and A = R'^W'^. By apphcation of this one-loop BT relation, 
supplemented by integration-by-parts identities, the power of the polynomial V(x) of the 
Feynman integral is raised. For example, for the one-loop three-point function in three 
dimensions, one finds 

rl pi— XI /'1—X1—X2 

/ dxi / dx2 / dxs Q(x)V'^(x) 

Jo Jo Jo 



pi pl—Xl pl—X\—X2 1 

- / dXi dX2 / dX3-r-V'^+^(x) 

JQ Jo Jo A 



' d 



2(/x + 1) ^ axfe 



fe=l 



/■I /■1-X2 1 1 

+ ^'A2arTT)-^'«x)V-(x)|_.. (12) 

This step can be applied iteratively until the power of V is as high as required, optimally 
II — —n+e — > e, where n is a positive integer. In this way, the original integral is expressed 
in terms of a sum of different integrals, which possess better integration properties. 

Although no general and compact-form solution of the BT algorithm is known for 
problems beyond the one-loop case, it is only natural to apply a one-loop BT relation to a 
sub-loop of a two-loop integral. In this way, the integral can be made smooth with respect 
to the Feynman parameters connected with the sub-loop to which the BT procedure is 
applied. Due to the size of the expressions and their divergency structure, it is usually 
better to apply this relation to the sub-loop with the highest number of internal lines. 

Finally, the e poles in each component Xi of the vector x are extracted by the relation 

'dx,x;-+j{^) 



I 



'0 

where i — 1, . . . ,n. 



k=0 / k=0 ^ ' 
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2.4 Semi-numerical integration based on dispersion relations 



This method makes use of the fact that the one-loop self-energy can be written, with the 
help of a dispersion relation, as an integral over an expression that has the analytical form 
of a propagator 



POO 

J (m-i +mo')2 



,2 



{mi+m2)^ ^ P 

where A(x, y, z) = + y"^ + — 2{xy + yz + zx). If this one-loop self-energy is a sub-loop 
of a two-loop integral, the dispersion relation effectively reduces this integral to a one-loop 
integral with the additional integration over s, which is performed numerically [21]. 

Integrals with sub-loop triangles can be reduced to integrals with sub-loop self-energies 
by introducing Feynman parameters [26]. The integration over the Feynman parameters 
is also performed numerically. More details can be found in Ref. [7]. If the two-loop 
integrals contains ultraviolet, infrared or threshold divergences, they need to be subtracted 
before the numerical integration can be carried out. For our purposes, this is achieved by 
subtracting a term from the integrand that can be integrated analytically. 

The reduction of integrals with irreducible numerators to a small set of master inte- 
grals is accomplished by using integration-by-parts and Lorentz-invariance identities. For 
complex diagrams with triangle sub-loops, the number of required relations can become 
very large, which is a limitation of this method. Therefore we do not use it to compute 
the complete result for the two-loop corrections to sin^ but only for cross checks of 
individual integrals and diagrams. 



2.5 Comparison of methods 

In this section, we compare our three methods, based on the top-quark mass expansion 
algorithm, the BT method, and the dispersion relations, where applicable. For the com- 
parison, we use the following dimensionless input parameters: Mz = 1, Mw = 80/91, 
rrit = 180/91, and the scale for dimensional regularisation fi = e"''^. The large- mass ex- 
pansion is performed down to (9(m^^^). The expressions are normalised as they enter 
in sin^6'gg, with the common prefactor (a;/47r)^ factored out. For the comparison, we 
selected the set of diagrams corresponding to Fig. [H where only IV-boson propagators in 
Feynman gauge, but not Goldstone-boson exchange has been included. For light-fermion 
loops, results are shown for leptons in the sub-loops, summed over the three lepton fami- 
lies, and for diagram (G) we chose the llv sub-loop, also summed over generations. Our 
results are presented in Table [H 

Where available, the results from the BT method and the method based on dispersion 
relations agree to all digit shown in the table. As mentioned above, no complete results 
for diagrams (E)-(G) could be obtained with the dispersion relation method. 
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Diagram Method Result [(a/47r)^] 

JA) semi-numerical +3.82775120/e^ +3.88128795/e -19.1983330 

mt expansion +3.82775120/e^ +3.88/e -19.19 

■(B) semi-numerical +3.82775120/e^ -8.67823832/e +25.4468576 

mt expansion +3.82775120/e^ -8.68/e +25.45 

JC) semi-numerical +0.90521614/e -0.60580110 

rrit expansion 0/e^ +0.905/e —0.61 

■(D) semi-numerical 0/e^ +1.55085212/e -4.50488822 

rrit expansion 0/e^ +1.55/e —4.50 

■(E) semi-numerical -2.30183413/e^ +5.07108758/e +8.32594367 

mt expansion -2.30183413/e^ +5.07/e +8.33 

■(F) semi-numerical -2.80183413/e^ +6.17261951/e -14.028 

mt expansion -2.80183413/6^ +6.17/e -14.03 

■(G) semi-numerical -1.80183413/e^ +3.9695556(5)/e -13.539 

mt expansion -1.80183413/e2 +3.97/e -13.54 



Table 1: Comparison of the top-quark mass expansion with semi-numerical integrations 
for selected diagrams of Fig. [H For diagrams (A)-(D), the semi-numerical results for the 
BT and dispersion relation methods agree to all given digits, while for diagrams (E)-(G) 
results are available only for the BT method. 

There are clear advantages to the use of large-m^ asymptotic expansions. No special 
treatment is required for the different types of divergences. Each large-mass pattern is 
associated with one expansion scheme, which is not sensitive to threshold divergences 
produced by the presence of small masses. In effect, the resulting programs are general 
and concise. As can be seen from Table [T], the obvious drawback of this method is the 
limited precision of the final results. However, it was observed that, with reasonable 
investment of computer time, the precision can be pushed to as high as required by the 
problem at hand. 

The semi-numerical programs usually produce results of better precision. However, 
the iterative application of the BT relation should be kept at a minimum, as otherwise the 
precision can actually be lost. At least in our realization, the semi-numerical programs 
based on the BT and dispersion relation methods are not as general as the expansion 
technique. Special care has to be taken to deal with threshold-divergent cases, and the 
size of intermediate expressions can actually be much larger than what is observed dur- 
ing the use of the large-mass expansion. In the end, the time required for the analytic 
simplifications of these semi-numerical programs can be as large as the time required for 
performing a high-order large-mass expansion. 

It should also be kept in mind that the algorithms for asymptotic expansions can be 
easily generalised to higher orders in perturbation theory, and were already applied for 
three-loop problems; see, for example, Refs. [27,28]. The application of the BT method 
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input pdl dlllfj Lt3I 


Value 




(HO 4-04 + n nnsn'i OpV 




(91.1876 ±0.0021) GeV 




2.4952 GeV 




(172.5 ±2.3) GeV 


rrib 


4.85 GeV 


Aa{Mz) 


0.05907 ±0.00036 


asiMz) 


0.119 ±0.002 




1.16637 X 10-5 GeV"2 



Table 2: Experimental input parameters used in the numerical evaluation, from Refs. [31, 
32]. 

in such cases is more complicated, and no physical results are known at this moment. 

3 Results 

The computational methods described in the previous section were implemented in au- 
tomatised computer codes, to be able to handle the large-size expressions at intermediate 
stages. The analytic algorithm for the BT reduction was written in FORM [29] and 
Mathematica. The code for the large-top-quark-mass expansion was also implemented 
in FORM. When necessary, the reduction of two-loop propagators with higher powers of 
numerators and denominators, which inevitably appear for higher orders of the asymp- 
totic expansion, was performed with the program IdSolver [4,25]. For the problem at 
hand, it had to create and solve a set of several thousands of equations, which it can 
achieve with very little computing time. For the numerical integrations, we developed a 
fast code written in C with the help of the Cuba library [30]. 

The results presented here were tested at different levels. We checked the finiteness 
and gauge invariance of the two-loop contributions to sin^ O^'^ analytically. The numerical 
results for the new diagrams were tested with two different methods, as presented in the 
previous section. In addition, full evaluation was performed independently by two groups 
within our collaboration. 

In the following, we show our numerical results for the input parameters listed in 
Table [2l For the sake of easy comparison, we use the same parameters as in previous 
publications on sin"^ 6^^^^ [7], which is justified by the fact that the changes of measured 
values are insignificant for this presentation. It is important to note that the experimental 
values for the W- and Z-boson masses correspond to a Breit-Wigner parametrisation with 
a running width and have to be translated to the pole-mass scheme used in the loop 
calculations [14]. In effect, this translation results in a downward shift [33] of Mz and 
Mw by about 34 and 28 MeV, respectively. The non-zero mass of the bottom quark 
was retained in the 0{a) contribution, but neglected in the two-loop part. The complete 
fermionic two-loop contribution to sin^ 9^^ is shown in Fig. [2], for various levels of the 
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Figure 2: Two- loop fermionic contribution to sin^ 9^^, with the W^-boson mass evaluated 
from the Standard Model and the other input parameters taken from Table El 



[GeV] [10-^] [IQ-^I [GeV] [10'^] [IQ-^J 



100 104.77 


1.00 


100 


105.03 


0.98 


200 100.15 


0.71 


200 


100.74 


0.67 


400 94.397 


0.32 


400 


95.354 


0.24 


600 90.666 


0.19 


600 


91.847 


0.10 


1000 85.748 


0.26 


1000 


87.196 


0.16 



Table 3: Contributions to sin^6'gg, with the VT-boson mass evaluated from the Standard 
Model (left) or fixed (right). 

large-top-quark-mass expansion. For the plot, is used as an input parameter, and 
Mw is calculated from it according to Eq. Q. As can be seen from the figure, the 
asymptotic expansion is converging robustly over the entire range of relevant Higgs-boson 
mass values. The relative error estimated for the expansion up to 0(m7^°) is 10-5 and 
thus more than sufficient for our purposes. The numerical values for selected values of the 
Higgs-boson mass are shown in Table [3l The left table uses as an input via Eq. (I3l), 
while a fixed mass Mw = 80.404 GeV is used for the right table. For small Higgs-boson 
mass, the new correction is relatively large, about 10"^. For larger values of Mh, it can 
be up to one order of magnitude smaller, about 10"^. Following earlier publications on 
two-loop electroweak corrections, we express our results in terms of fitting formulas. The 
form factor Ak^^ ,ferm) ^j^j^,]^ contains the fermionic two- loop electroweak corrections to 
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sin according to Eq. ([2]), can be approximated as 

00 00 00, rem 

^4£em"^ = ^0 + ^^^^ + ^2^^ + ^3^^ + ^^l^H " 1) + 



+ hAf + kjAtLn + ksAw + kgAwAt + ki^Az, (15) 



where Ak"^ is the one-loop result, and 

^^-^^lOOGiV' ^^-lOOG^' - llTSGivJ " ^' 

- 91.1876 GeV " ^' " 80.404 GeV ~ ^^^^ 

Fitting this formula to the exact result, we obtain 

ko = -2.666 X 10"^ ki = -5.92 x 10~^ ^2 = -3.29 x 10"^ 

fcg = 3.49 X 10"^ A;4 = 2.83 x 10~^ fcg = -5.34 x 10"^ 



fee = -2.10 X 10"^ A;7 = -2.19 X 10^^ fcg = -6.31 x 10"^ 
fcg = -1.26 X 10~\ kw = 6.47 X 10"=^. 



(17) 



This parametrisation reproduces the exact calculation with maximal and average devia- 
tions of 1.4 X 10~^ and 5 x 10~^, respectively, as long as the input parameters stay within 
their 2a ranges of the experimental errors quoted in Table [2] and the Higgs-boson mass is 
in the range 10 GeV < Mh < 1 TeV. If the top-quark mass and the H^-boson mass vary 
within 4(7 ranges, the formula is still accurate to 2.1 x 10~^. 

We also present a simple parametrisation for the currently best prediction for sin^ O^g, 
including all known corrections to Ak^^ and Ar (for the calculation of from see 
Refs. [2,34]). For Ak^^, in addition to the one-loop and fermionic two-loop electroweak 
corrections, we include QCD corrections of 0{aas) [35] and O(aa^) [27,36] to the one- 
loop contribution, as well as universal corrections for large top-quark mass, of 0{a'^asmf) 
and 0{a^m^) [28]. Moreover, leading four-loop QCD correction to the p parameter, 
which arise from top- and bottom-quark loops, are taken into account [37]. We use the 
parametrisation 

sin^ = So + diLu + ^2^^ + + 4(A^ - 1) + c^sA^ 

+ deAt + djA^t + dsAt{AH - 1) + dgA^, + dioAz, (18) 

with 

0.05907 ' 0.117 ^ ^ 

The best-fit numerical values for the coefficients are 

So = 2.327580 x 10-\ di = 4.749 x 10'^, d2 = 2.03 x 10-^ 



4 = 3.94x10"^ 4 = -1.84 X 10"^ rfg = 2.08x^°"2, 

4 = -9.93 X 10-^ dj = 7.08 x 10"^ 4 = -7.61 x 10" 
dg = 4.03 X 10-^ dio = 6.61 x 10"\ 



(20) 
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This parametrisation approximates the full result with maximal and average deviations 
of 4.3 X 10"^ and 1.3 x 10"^ respectively, for 10 GeV < Mh < I TeV and the other input 
parameters in their 2a ranges. 

4 Conclusions 

In this paper, the calculation of the two-loop electroweak fermionic corrections to the 
effective weak-mixing angle for the Zbb vertex, sin^ 9^^, was presented. Such an accurate 
theoretical prediction for sin^ 9^^ is necessary for the interpretation of the bottom-quark 
asymmetry measurements at the Z-boson pole. Compared to the previously known cor- 
rections to sin^^gg, the new electroweak two-loop result turns out to be sizable, of order 
(9(10~^) for a Higgs-boson mass near 100 GeV. 

The calculation was performed by using methods that had been used earlier for the 
computation of the leptonic effective weak-mixing angle, as well as a newly developed code 
based on the BT algorithm. The results of the different methods were checked against 
each other. 

Although we did not perform a detailed analysis of the error from unknown high-order 
corrections, in particular the missing bosonic two-loop corrections and terms of order 
0{a'^as), we expect those to be of similar order as for the leptonic effective weak-mixing 
angle. The main difference between the leptonic and bottom-quark effective weak-mixing 
angles are the vertex diagrams with internal VT-boson and top-quark propagators. While 
leading to numerical differences between sin^ 6'''^* and sin^^^l, these diagrams do not 
introduce special enhancement or suppression factors. Therefore, we expect the theoretical 
uncertainty to our result for sin^ 9^^ to be about 5 x 10"^, similar to Ref. [7]. 
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